2HDM Portal Dark Matter: LHC data and the Fermi-LAT 135 GeV Line 
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We study a two Higgs doublet model augmented by a scalar dark matter particle that provides 
an excellent fit to the LHC Higgs data and the Fermi-LAT 135 GeV line. The heavy CP-even Higgs 
boson, which predominantly mediates annihilation and scattering, must have a coupling to weak 
gauge bosons at or below percent level to suppress the continuum gamma-ray spectrum below the 
limit from the Fermi-LAT data and the anti-proton spectrum constrained by the PAMELA data. 
Discovering or excluding this CP-even Higgs boson at the LHC with a mass between 265 and 280 
GeV and an enhanced diphoton branching ratio is crucial to test this scenario. 
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Introduction. The way in which dark matter interacts 
with the Standard Model (SM) remains a mystery. Sce- 
narios in which the photon or Z boson are mediators with 
electroweak interaction strengths have already been ex- 
cluded by dark matter direct detection experiments. The 
remaining natural mediator is the SM Higgs boson, re- 
sulting in the "Higgs portal" dark matter scenario [l|43| ■ 
With the recent discovery of a Higgs-like particle with a 
mass of 125-126 GeV at the LHC @, i], the SM Higgs 
portal dark matter scenario is excluded for dark matter 
masses below 1 TeV except for the resonant region in 
which the dark matter mass is close to one half of the 
Higgs mass @ . This suggests that extensions to the sim- 
plest Higgs portal are worthy of further consideration. 

The recent LHC Higgs data, which includes a rela- 
tively large diphoton branching ratio, also hints more 
generally at new physics in the Higgs sector. A simple 
way to extend the SM Higgs sector is to consider models 
with two Higgs doublets (2HDM's) 0-0. In such scenar- 
ios, the mixture of the two CP-even neutral scalar fields 
and general couplings of two Higgs doublets to fermions 
can dramatically modify the lightest Higgs boson prop- 
erties The increase of the light Higgs diphoton 
branching ratio in the generalized 2HDM, which we have 
called the 2HDM-X Q, can be realized in two ways: re- 
ducing the total width or increasing the coupling to two 
photons from the charged Higgs contributions. 

In addition, there has recently been a hint of a gamma- 
ray line around 130 GeV from the galactic center that 
was found from analyzing the Fermi Gamma-Ray Space 
Telescope (Fermi-LAT) data ^iMiM- This hint has 
been confirmed by the Fermi-LAT collaboration with a 
smaller statistical significance. Using the reprocessed 
data, the peak has shifted to a slightly higher mass at ~ 
135 GeV [16[. Many dark matter models with additional 
charged particles have been constructed to explain this 
tentative Fermi-LAT 135 GeV gamma-ray line (see [13] 
and references therein). 

In this letter, we propose a "2HDM portal" dark mat- 
ter scenario that can explain both the LHC Higgs data 



and the Fermi-LAT 135 GeV gamma-ray line. The char- 
acteristic feature of the scenario is that the heavy CP- 
even Higgs boson of the 2HDM-X model is the mediator 
for the SM particles interacting with the dark matter 
particles. For a wide range of parameter space, the dark 
matter particles primarily annihilate into two gamma-ray 
lines with suppressed cross sections into the continuous 
gamma-ray and anti- proton spectra [l8l - l20| . 

Motivated by the possibility of realizing the first order 
electroweak phase transition for baryogenesis [21 , 22 1 , we 
consider the dark matter to be a real SM singlet scalar 
field in this work, and defer the consideration of the case 
of fermionic dark matter for a future study. For the case 
of SM singlet scalar dark matter, the interactions of the 
dark matter can be described by a rcnormalizablc scalar 
potential with only a few new parameters. We will sec 
that for this scenario to accommodate both the collider 
and astrophysical data, the heavy neutral Higgs boson 
must have essentially a vanishing vector boson coupling, 
though with a strict mass range between 265-280 GeV. 

The 2HDM portal. We consider the two Higgs doublet 
model of [|| (the 2HDM-X), which includes two complex 
scalar doublets of opposite hypercharge: 
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(246 GeV) 2 , and the ratio 



of the two Higgs vacuum expectation values is tan /? = 

To extend the 2HDM-X to include dark matter, we 
include a real scalar field, S, which is assumed to obey a 
Z2 symmetry (S — s- — S). 1 In the absence of CP violation, 



1 One can also consider a complex scalar, where the CP-odd com- 
ponent can be an equally viable dark matter candidate l23l . How- 
ever, much of the phenomenological discussion can be mapped 
to the real scalar case, which we focus on here. 
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the scalar potential is given by 

V = V 9 + Vs + V 9S , (2) 
where V$ is the CP-conserving 2HDM scalar potential, 

V$ = m 2 <f>|$i + m%<b\<&% - [m\ 2 <&\<& 2 + h.c. 
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(*i,2 = -**,2 icr 2), Vs is given by 
V S 
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and V$s, which includes the interactions between the 
scalar singlet and the Higgs doublets, takes the form 

8 8 8 

V<S, S = yS 2 $t$l+yS 2 *£$ 2 + -!S 2 ($I$2+h.C.), ( 5 ) 

in which S3 is taken to be real. As we will see, the cou- 
plings #1,2,3 are responsible for many observable effects. 

The Yukawa interactions of the Higgs doublets with the 
SM fcrmions are assumed to take the following restricted 
form, as discussed in Q: 

- C = y u u R (cos 7„ * 2 ~ sin j u #1) Q L 
+ ydd R (cos7 d *i + sin7 rf * 2 )QL 
+ yi~e R *i L L + h.c. , (6) 

in which the y u ,d,i are 3x3 matrices in family space, and 
the j Ut d lie within the range < j U)d < n. 

At leading order, the CP-even ^77 effective couplings 
are given by 
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where in the sum over each loop particle j, N C j is the 
color factor, Qj is the charge, Kij is the coupling to hi 
relative to the SM coupling, the Fj are the respective 
loop functions, and = 4m| /m^ . . The gluon am plit ude 
follows similarly, but retains only the quark loop 24 1 . 

The tree-level deviations of the Yukawa couplings of 
the neutral Higgs bosons from the SM Higgs boson are 
given in Table U and will generally referred to as ftj, with 
mass index i. The CP-odd and charged Higgs couplings 
scale with A u = cot(/3 + 7„), Ad = tan(/3 — 7^) and 
A( = tan/3. The charged Higgs couplings are given by 

V2 
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TABLE I: The tree level Yukawa couplings of the CP-even 
(h, H) and CP-odd (A) states with respect to the SM Higgs. 
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We exchange the dimensionlcss scalar potential couplings 
for the masses (Mh,MH,MA,M H ±,Ms) and mixings 
(the usual CP-even Higgs mixing angle a and tan j3) 
while maintaining pcrturbativity of the scalar couplings. 

Dark Matter Observables. The direct detection of 
dark matter is observed via scattering off nuclei, which 
proceeds in this model through an exchange of neutral 
Higgs bosons, as shown in Fig. [TJ The proton- S spin- 




FIG. 1: Feynman diagrams for the scattering of dark matter 
off the partons of the proton. 



independent scattering rate is given by 
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in which the index i accounts for both CP-even Higgs 
states. The SShi couplings, which arise from V$s, take 
the following values: 
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- 82 cos a sin (3 — 8\ sin a cos /3 + 83 cos(a + (3) , 
= 81 cos a cos (i + 82 sin a sin /3 + 83 sin(a + /3)(11) 



Indirect detection begins with dark matter annihilation. 
Dark matter in the present epoch is non-relativistic, and 
hence we treat annihilation as occurring in the static limit 
v — > 0. Since here the dark matter is a scalar field that 
couples only to the Higgs boson, one can write a portion 
of dark matter annihilation rate in terms of the Higgs 
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provides the following distilled measurements [3. [HI. |2£ 
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FIG. 2: Feynman diagrams for the annihilation of dark matter 
to photon pairs. 



decay rate, as follows: 



hSM^XX 



x( m h s 



2m s ) 



ss^xx u — 
gssh ki(XX) 



ims 



(12) 



gSSH K2(XX) 



4m| - m 2 h 



im h T h 



4m| 



im H T H 



in which X denotes any SM state that is not a Higgs 
boson or photon. For photons, the annihilation rate is 
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which contains the processes shown in Fig. [21 The 
remaining annihilation modes include SS — > hh and 
hH, HH, AA, H + H~ when kinematically accessible. The 
Higgs bosons may acquire an additional decay width of 



T(hi -> SS) = 



9ssh z 
32iTmh i 



(14) 



which may be seen at colliders as missing energy Indeed, 
searches for invisibly decaying Higgs bosons may be 
possible with early LHC data [25| ■ 

Fit to available data. We fit the 2HDM portal sce- 
nario outlined above to the available data using Baycsian 
inference, specifically adopting the Metropolis-Hastings 
Markov Chain Monte-Carlo algorithm; details for our im- 
plementation of this algorithm are given in [26j . Among 
its advantages include an efficient exploration of the 
model space as well as readily available posterior proba- 
bility distributions. In addition to applying the experi- 
mental constraints listed below in our fit, we require all 
couplings remain perturbative: g,S, A < \/4tt. 
• Collider Data. We perform a Bayesian fit to the avail- 
able collider data following the same method as 27j . We 
require that the lightest CP-even mass eigenstate is the 
boson observed at the LHC and the Tevatron. A simple 
combination of inclusive and selected exclusive channels 



We also require the second CP-even Higgs state to fall 
below the present exclusion limits from ATLAS, which 
strongly disfavor a large vector boson coupling for moder- 
ate Higgs masses. In addition, we include limits from the 
flavor changing decays Bo — > X s +-f+X at next to leading 
order [3~fl ] for consistency with the experimental measure- 
ment of BF(S -> X s +7+ X) = (3.55±0.26)xl0" 4 [H]. 
We further require consistency of the electroweak oblique 
parameters with S = 0.04 ± 0.09, T = 0.07 ± 0.08 with 
an 88% positive correlation |33| . Recent measurements 
of BF(£? S -+ /i + /i~) have little impact on our scan as the 
new physics effects are suppressed for low tan/3. 
• Dark Matter Data. We assume the Fermi 77 line at 
E-y = 135 GeV is astrophysical in origin and arises from 



dark matter annihilation with ms 
the fitted value UM 



135 GeV. We adopt 



a^v = 2.27t°; 6 ^ (1.27t£S) X 10 



+0.37^ 



1-27 cm 3 s _1 , 



(20) 



which assumes an NFW (Einasto) dark matter galactic 
halo profile. The astrophysical uncertainties for the re- 
quired dark matter annihilation cross sections are fairly 
large. While we include the NFW value in our fit, we 
show the Einasto fit in our results as a comparison. Gen- 
erally, since the Einasto profile yields a lower fitted cross 
section, it can be more easily accommodated within this 
model. The Fermi collaboration analysis finds a simi- 
lar line, but with lower significance. In addition to the 
SS — > 77 line, a generic model will have annihilation to 
SS — > Z7, where the 7 line is shifted to lower energy. 
In this scenario, E 1 = 120 GeV. We impose a maximal 
cross section via Fermi data for this mode of II 611 



az-yV < 1.4 x 10 



-"emV 1 , 



(21) 



The Fermi satellite further includes limits on the sec- 
ondary photons from possible dark matter annihilation in 
dwarf galaxies and the galactic halo [2(|. Furthermore, 
stringent limits on annihilation to hadronic final states 
can be placed from PAMELA p data. For the mass we 
assume, the combined limits are strongest for annihila- 
tion to bb and W + W~ . We require the annihilation cross 
section for the most constrained modes to fall below 134] 



26 cm 3 s _1 , 



a bi v < 4.2 x 10 
cr w +w-v < 3.8 x 10 _26 cm 3 s- 1 
(? t +t-v < 1-4 x 10" 25 



(22) 
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Due to the dark matter mass indicated by the Fcrmi- 
LAT line, another important annihilation mode includes 
SS — > hh annihilation which can proceed via s-channel 
h/H, i-channel S and four-point SShh processes. This 
mode is taken into account by scaling the energy from 
dark matter annihilation by 1/2 to account for halv- 
ing the available phase space for the annihilation prod- 
ucts, while the cross section limit increases by a factor 
weighted by the respective branching fraction for the final 
state particles. For instance, for bb, we have 



<Thh{m s ) - v = <7 b - b ■ v 



1 



2BF(/i^ bb) 



(23) 



where either or both h's can decay to bb. We also apply 
the most restrictive direct detection measurement to 
date for dark matter of moderate mass, which comes 
from the XENON experiment. This gives a limit of 
cr SI < 4 x 1(T 45 cm 2 for m s = 135 GeV at 95% C.L. 

Discussion. The LHC and Tevatron Higgs data ap- 
pear to prefer a hVV coupling consistent with custo- 
dial symmetry, since both the ratio in the WW/ZZ and 
overall rates agree within measurement uncertainty with 
the SM. These constraints, coupled with the null re- 
sults for heavier Higgs masses, leave little room for a 
large HVV coupling. Therefore, the key observation is 
that if annihilation is mediated primarily through H , the 
H — > W + W~~ , ZZ modes are suppressed, allowing for an 
escape from the limits in Eq. (|22j) . Such a dominant an- 
nihilation through H occurs when the SSh coupling is 
suppressed, which happens when 
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where Ay, which is defined by g\ vv = 1 
tcrizes the level of decoupling in the heavy Higgs sector. 

We also note that a suppression in the HVV coupling 
from the SM allows a natural enhancement of the SS — > 
77 line relative to W + W~ . Indeed, the 77 line fit and 
continuum 7 ray constraints require at least 



BF(</> -> 77) 



BF( 



w+w-) 



— > e>(l(T 2 ) ; 



(25) 



for = 2ms = 270 GeV. In contrast, a SM Higgs 
mediating dark matter annihilation gives only 



BF(0 77) 
BF(0^ W+W') 



0(icr 5 ). 



(26) 



This enhancement thus does not occur in a simpler Higgs 
portal model with only one Higgs boson. 

The SS — ¥ hh annihilation mode is kincmatically ac- 
cessible and is of importance due to the h decay to sec- 



ondary 7 and p. 2 The constraints on this mode require 
a suppression in both gsshh and gssh couplings, which 
take similar form gssh ~ v gsshh- Therefore, the sup- 
pression required for the SS — > bb and W + W~ channels 
also partially suppresses the SS — > hh amplitudes in- 
volving the s-channel h and four-point SShh diagrams. 
The remaining s-channel H diagram is suppressed in the 
decoupling limit by the coupling 



ghhH — Ai 
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FIG. 3: The correlation of the decoupling limit parameter 
Ay with a 11 v, which shows the importance of the suppressed 
WWH and ZZH couplings in fitting the Fermi-LAT data. 
The vertical lines illustrate the la range for the Fermi-LAT 
line. 



The decoupling limit and the constraints from sec- 
ondary 7 and p thus provide an explanation for the size 
of the Fermi-LAT 77 line within this scenario. The im- 
portance of the decoupling limit is highlighted in Fig. [31 
which shows that to reproduce the Fermi-LAT line, the 
decoupling limit must be realized to better than the mil- 
level. 

Due to the tight connection between the heavy Higgs 
boson and the associated annihilation cross sections, the 
Fermi-LAT line may also provide insight to the heavy 
Higgs mass. In Fig. HI we observe a tight region for run 
that centers on the resonance SS — > H from 265-280 
GeV at 95% C.L., but with a narrow gap at 270 GeV. In 



2 We enforce similar restrictions to the SS — ¥ H + H~ and Hh 
modes, when accessible. The light charged Higgs is preferred to 
give an enhanced h — > 77 rate seen at the LHC. The AA mode 
can easily be closed by increasing rax above threshold. 
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FIG. 4: The correlation of the heavy and charged Higgs 
masses with a 11 v. To explain the Fermi-LAT line, the heavy 
Higgs mass is likely to be near the SS — > H resonance. The 
vertical lines illustrate the la range for the Fermi-LAT line. 



this region, the ii-rcsonant annihilation to bb increases 
beyond the Fermi limit. The charged Higgs mass is 
expected to lie in a region of 140 — 320 GeV at 95% 
C.L., as required to elevate the h — > 77 and SS — > 77 
rates. We also note that if the 77 rate seen at the LHC 
subsides to the SM-like rate, the expected range for run 
and tojj± does not appreciably change. 

Summary. The 2HDM portal dark matter scenario, in 
which the 2HDM is extended with a real scalar singlet 
that plays the role of the dark matter, can explain both 
the LHC Higgs data and the Fermi-LAT gamma ray line 
data. In this model, the heavy CP-even Higgs state is 
the mediator of the interactions of the dark matter with 
the SM fields. The model predicts the following Higgs 



sector features: (i) the heavy Higgs sector is decoupled 
at the percent level or better, (ii) the heavy CP-even 
Higgs mass is preferentially between 265-280 GeV, and 
(iii) the charged Higgs mass is preferentially between 140 
and 320 GeV. Furthermore, the production of secondary 
7 and p are suppressed below present bounds, but it is 
not unrealistic to expect an observation of these modes 
soon as the SS — > hh annihilation rate can be large. 
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